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Abstract 

We consider continuous state branching processes that are perturbed by a Brownian 
motion. These processes are constructed as the unique strong solution of a stochastic 
differential equation. The long-term extinction and explosion behaviours are studied. 
In the stable case, the extinction and explosion probabilities are given explicitly. We 
find three regimes for the asymptotic behaviour of the explosion probability and, as 
in the case of branching processes in random environment, we find five regimes for 
the asymptotic behaviour of the extinction probability. In the supercritical regime, 
we study the process conditioned on eventual extinction where three regimes for 
the asymptotic behaviour of the extinction probability appear. Finally, the process 
conditioned on non-extinction and the process with immigration are given. 
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1 Introduction 

A [0, oo]-valued strong Markov process Y = (Vt,t > 0) with probabilities (Pa,, x > 0) is 
called a continuous-state branching process (CB-processes for short) if it has paths that 
are right-continuous with left limits and its law observes the branching property; i.e. for 
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any x, ?/ > 0, is equal in law to the convolution of P^; and P^. CB-processes may 
be thought of as the continuous (in time and space) analogues of classical Bienayme- 
Galton-Watson processes. CB-processes have been introduced by Jirina [13] and studied 
by many authors including Bingham [5], Grey m, Lamperti [18] , to name but a few. The 
branching property implies that the Laplace transform of satishes 


Ex 





for A > 0, 


( 1 ) 


for some function Ut{X). According to Silverstein [23], the function Ut{\) : [0, oo) —?■ [0, oo] 
solves the integral equation 

Ut{X) + [ 'ilj{us{X))ds = X, ( 2 ) 

Jo 

where ^ satishes the celebrated Levy-Khintchine formula, i.e. 

'ijj{X) =-q - aX ++ f - 1-h Aa;l{ 3 ;<i})/i(da;), 

J (0,oo) 

where g>0, ogM, 7>0 and /r is a measure concentrated on (0, cxo) such that (l A 

x^)qi{dx) is hnite. The function ijj : [0, cxo) {—oo,oo) is convex and is known as the 
branching mechanism of Y. 

Observe that 0 and oo are two absorbing states. In other words, let 


To = inf{f > 0 : Tt = 0} and Too = inf{t > 0 : Yf = oo} 


denote the extinction and explosion times, respectively. Then 1} = 0 for every t > Tq 
and Yt = oo for every t > Too- More precisely, let rj be the largest root of the branching 
mechanism ip, i.e. 7] = sup{0 > 0 : ip{9) = 0}, (with the convention that sup{0} = oo). 
Then for every x > 0: 


i) if ?7 = oo or if d9/ip{9) = oo, we have Px(To < oo) = 0, 

ii) if rj < oo and /°°d6'/'0(0) < oo, we dehne 


0(f) = 


d9 

wr 


t e (r/,oo). 


The mapping 0 : {rj, oo) —?■ (0, oo) is bijective, and we write (p : (0, oo) — {rj, oo) for 
its right-continuous inverse. Thus 


P^(To <t) = exp{-x</?(f)}. 


hi) if 7 = 0 or if d6'/|0(6*)| = oo, we have Px(Too < oo) = 0, 
iv) if ?7 > 0 and /g_|_ d6'/|■0(0)1 < oo, we dehne 

. ^ f* d9 , 
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The mapping g : {0,r]) —)■ ( 0 , cxo) is bijective, we write 7 : ( 0 , cxo) ( 0 , 7 ) for its 
right-continuous inverse. Thus 

^xiToo > t) = exp{-a:7(t)}. 

From (ii), we deduce that Px(To < 00 ) = exp{—Hence, the latter identity and (i) 
imply that a CB-process has a hnite time extinction a.s. if and only if 

du 

n < 00 , / , . . < 00 and '0'(O-|-) > 0. 

Similarly from (iv), we get that P 3 ;(Too < cxo) = 1 — exp{—xr^}. Hence from the latter and 
(hi), we deduce that a CB-process has a hnite time explosion with positive probability if 
and only if 

f 

/ I <00 and 7 > 0, 

Jo+ mu)\ 

When ?7 < 00 , 7 > 0 is equivalent to < 0. 

The value of -^'(O-l-) also determines whether its associated CB-process will, on average, 
decrease, remain constant or increase. More precisely, under the assumption that q = 0, 
we observe that the hrst moment of a CB-process can be obtained by differentiating ([T]) 
with respect to A. In particular, we may deduce 

for t > 0 . 

Hence using the same terminology as for Bienayme-Galton-Watson processes, in respective 
order, a CB-process is called supercritical, critical or suhcritical depending on the behaviour 
of its mean, in other words on whether < 0 , = 0 or > 0 . 

The following two examples are of special interest in this paper since the Laplace 
exponent, i.e. the solution of ([2]), can be computed explicitly in a closed form. The 
hrst example that we present is the so-called Neveu branching process (see [SU]) whose 
branching mechanism satishes 

f cIt* 

V'(A) = Alog(A) = cA- 1-F Axl{a;<i})—, 

J ( 0 ,oo) ^ 

where c G M is a suitable constant. In this case -^'(O-l-) = — 00 , 7 = 1 , the process 
is supercritical and satishes the integral conditions of (i) and (iii). Thus, the Neveu 
branching process does not explode neither become extinct at a hnite time a.s. According 
to Theorem 12.7 in na, we have that the Neveu branching process becomes extinct at 
inhnity with positive probability. More precisely, if (T), f > 0) denotes the Neveu branching 
process, then 

Pa; Qim Ft = 0j = e~^, x > 0. 

The second example is the stable case with drift, in other words the branching mechanism 
satishes 

= —a\+cpX^^^ = —a\+C2\^l{li=i}+c^[^ - f (e l-|-Axl{y3>o}) 2+^ 

tl-L-pj J( 0 ,oo) ^ ^ 
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with jS in (—1, 0) U (0,1] and is a non-zero constant with the same sign as 13. It is known 
that its associated CB-process can be obtained by scaling limits of Bienayme-Galton- 
Watson processes with a hxed reproduction law. Moreover, the case (3 = 1 corresponds to 
the so-called Feller diffusion branching process. 

The case (3 G (—1,0) has a particular behaviour. Its corresponding CB-process is 
supercritical, since V^'(O-I-) = —oo. We observe that rj is inhnite or hnite according to 
whether a > 0 or a < 0. The process satishes the integral conditions of (i) and (iv). Thus 
the stable CB process with (3 G (—1, 0) does not become extinct at a hnite time a.s., and it 
has a hnite time explosion with positive probability. Moreover, the asymptotic behaviour 
of P 3 ;(T'oo > t) can be computed explicitly. 

The case [3 G (0,1] has a completely diherent behaviour. Its associated CB-process 
is subcritical, critical or supercritical depending on the value of a, since '0'(O-|-) = —a, 
and satishes the integral conditions in (ii) and (iii). We also have that 77 = 0 or 77 > 0 
according as a < 0 or a > 0. In other words, the stable CB-process with (3 G (0,1] does 
not explode at a hnite time a.s., and it becomes extinct at a hnite time with positive 
probability. Moreover the asymptotic behaviour of Pa;(To < t) can be computed explicitly. 

For our purposes, we recall that CB-processes can also be dehned as the unique non¬ 
negative strong solution of the following stochastic diherential equation (SDE for short) 


Yt =Yo + a Wds 



r-y,- 


2 ;A^(ds, dz, du) -|- 


0 ^( 0 , 1 ) Jo 



0 ./[l,oo] Jo 


-Ys- (3) 

zA^(ds, dz, du). 


where B = {B^, t > 0) is a standard Brownian motion, N{ds, dz, du) is a Poisson random 
measure independent of B, with intensity dsA(dz)dM where A is a measure on (0, 00 ] de¬ 
hned as A(d 2 ) = l((o,oo)}(^)Ai(cl^) + qSoo{dz), and N is the compensated measure of N, see 
for instance HDl. 


In this work, we are interested in studying a particular class of CB-processes in a 
random environment. More precisely, we are interested in the case when the random 
environment is driven by a Brownian motion which is independent of the dynamics of 
the original process. A process in this class is dehned as the unique non-negative strong 
solution of a stochastic diherential equation that conditioned on the environment, satishes 
the branching property. We will refer to such class of processes as CB-processes in a 
Brownian random environment. 

Our motivation comes from the work of Boinghoh and Hutzenthaler [ 6 ] and Hutzen- 
thaler [12], where they consider the case of branching dihusions in a Brownian random 
environment i.e. when the branching mechanism has no jump structure. The authors in 
[ 6 ] introduced this type of processes using a result of Kurtz [H], where a dihusion approxi¬ 
mation of BienaymCCalton-Watson in random environment is studied. The scaling limit 
obtained in [T3| turns out to be the strong solution of an SDE, that conditioned on the 
environment, satishes the branching property. Boinghoh and Hutzenthaler computed the 
exact asymptotic behaviour of the survival probability using a time change method and in 
consequence, they also described the so called Q-process. This is the process conditioned 
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to be never extinct. Similarly to the discrete case, the anthors in p] fonnd a phase transi¬ 
tion in the snbcritical regime that depends on the parameters of the random environment. 
Hntzenthaler stndied in [12], the snpercritical regime and fonnd that snpercritical bran¬ 
ching diffnsions in a Brownian random environment conditioned on eventnal extinction 
also possesses a phase transition which is similar to the phase transition of the snbcritical 
case. 

Another class of CB-processes in random environment has been stndied recently by 
Bansaye et ah [1]. The anthors in [1] stndied the particular case where the random 
environment is driven by a Levy process with paths of bounded variation. Such type 
of processes are called CB-processes with catastrophes, motivated by the fact that the 
presence of a negative jump in the random environment represents that a proportion of a 
population, following the dynamics of the CB-process, is killed. Similarly to the diffusion 
case, CB-processes with catastrophes were also introduced as the strong solution of an 
SDE and conditioned on the environment, they satisfy the branching property. It is also 
important to note that CB-processes with catastrophes can also be obtained as the scaling 
limit of BienaymAGalton-Watson processes in random environment (see for instance 0). 
Bansaye et ah also studied the survival probability but unlike the case studied in [B], they 
used a martingale technique since the time change technique does not hold in general. In 
the particular case where the branching mechanism is stable, the authors in [1] computed 
the exact asymptotic behaviour of the survival probability and obtained similar results to 
those found in [Bj. 

One of our aims is to study explosion and extinction probabilities for CB-processes in 
a Brownian random environment. Up to our knowledge, the explosion case has never been 
studied before even in the discrete setting. In order to study explosion and extinction 
probabilities for a process in this class, we follow the martingale technique used in [T] to 
compute the Laplace exponent via a backward differential equation. With the Laplace ex¬ 
ponent in hand, we are able to determine whether a process is conservative, i.e. that does 
not explode a.s. at a hxed time, or become extinct with positive probability. Nonethe¬ 
less, it seems very difficult to deduce necessary and sufficient conditions for explosion and 
extinction probabilities. This is due to the fluctuations of the random environment. How¬ 
ever, the stable and the Neveu cases will help us to understand the different situations 
that our main results cannot cover. 

We give special attention to the case when the branching mechanism is stable i.e. 
■^(A) = —aX + CfiX^^^, for jS G (—1, 0)U(0,1]. Here, the Laplace exponent can be computed 
explicitly and we will show that it depends on the exponential functional of the random 
environment. Whenever f3 G (—1,0), we can compute the explosion probability at a hxed 
time and establish the asymptotic behaviour of the probability of no-extinction where we 
hnd three different regimes. Up to our knowledge, this behaviour was never observed. 
In the case when /3 G (0,1], we study the extinction probability and also establish the 
asymptotic behaviour of the survival probability where hve different regimes appear. In the 
supercritical regime, we study the process conditioned on eventual extinction and we hnd 
three regimes for the asymptotic behaviour of the probability of survival. The asymptotic 
behaviour depends on the study of exponential functionals of Brownian motion. 

From the speed of survival in the stable case, we can deduce the process conditioned 
to be never extinct or Q-process using a Doob h-transform technique. 


5 


We finish this paper studying the immigration case, which represents an example of 
affine processes in a random environment. More precisely, this family of processes is an 
extension of the so-called Cox-Ingersoll-Ross model, which is largely used in the hnancial 
literature, under the fluctuations of a random environment. 

The remainder of the paper is structured as follows. In Section 2, we dehne and study 
CB-processes in a Brownian random environment. Section 3 is devoted to the long-term 
behaviour of this family of processes. In particular, we study explosion and extinction 
probabilities. In Section 4, we analyse the stable case. Here, we study the asymptotic 
behaviour of the no-explosion and survival probabilities as well as the process conditioned 
to be never extinct and the process conditioned on eventual extinction. Finally in Section 
5, we study the immigration case. 


2 CB-processes in a Brownian random environment 


Motivated by the dehnition of branching diffusions in random environment (see Boinghoff 
and Hutzenthaler 0 ) and CB-processes with catastrophes (see Bansaye et al. [1]) we 
introduce, using the same notation as in the SDE ([3]), continuous state branching processes 
in a Brownian random environment (in short a CBBRE) as the unique non-negative strong 
solution of the following stochastic differential equation 


Zi =Zo + a Z,ds + / V'V^dR, + a Z,dR(^) 



-Zs- 


zN{ds, dz, du) -|- 


0 ^(0,1) Jo 



(4) 


2 ;V(ds, dz, du). 


0 ^[l,oo] Jo 


where > 0) is a standard Brownian motion independent of B and the 

Poisson random measure N, and a G M. The Brownian motion Rh) represents the random 
environment. We also observe that the drift coefficient can be written as a = a + do, where 
a is the drift of the underlying CB-process and Oq is the drift of the environment. 

The following theorem provides the existence of the CBBRE as a strong solution of 
(l4l) and, in some sense, characterizes its law given the environment. In order to introduce 
our main result, we dehne the auxiliary process 

2 

Kt = (tB[^^ — —t, for t > 0, 

that represents the random environment. 

Theorem 1. The stochastic differential equation (j4]) has a unique non-negative strong 
solution. The process Z = {Zt,t > 0) satisfies the Markov property and its infinitesimal 
generator C satisfies, for every f G C^(M+)j3 


Cf(x) =axf{x) 


2 

-I- X 


2 2 I 2 

a X j X 


(0,oo] 


{f{x + z) - f{x) - zf{x)l{,<i}) A(d 2 r). 


( 5 ) 


t_|_ = [0,oo), IR+ = [0, oo] and (]R+)={twice differentiable functions such that /(oo) = 0} 
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Furthermore, the process Z conditioned on K, satisfies the branching property and for 
every t > 0 and z > 0 


E, 


exp 


XZtC 


-Kt 


K 


= exp 




a.s.. 


( 6 ) 


where for every (A, 5) G (M+,C'(M+)), Vt : s E [0,t] Vt{s,X,6) is the unique solution of 

the backward differential equation 

-^vt{s,X,6) = e^‘^ilj{vt{s,X,6)e~^‘’), vt{t,X,6) = X, (7) 

and if is the branching mechanism of the underlying CB-process. 


Proof. First, we prove the existence of a unique strong solution of the SDE 







An)d5,+ 



+ / / / (z A n)A^(ds, dz, du) + / / / A n)A^(ds, dz, du). 

Jo J{0,1) Jo Jo J[l,oo] Jo 

( 8 ) 


Let E = {1,2} and 7i{dz) = 5i(d^) + S 2 {dz) a measure in E. Dehne hF(ds,d^) := 
di?s(5i(dz) + dBi^'* 62 {dz), a white noise on [0, oo) x E with intensity ds7r(d2;). Then, the 
SDE (IH]) can be written as follows 

pin) f ^ 

Jo 

+ 11^ An) + A lE(ds, d^) 

rt f- i‘{zi”fin)- _ 7* 7 Azf^fin)- 

+ / (z A n)iV(ds, dz, du) + / / / (z A n)A^(ds, dz, du). 

Jo J{0,1) Jo Jo j[l,oo) Jo 

(9) 

Following the notation in [5], the conditions from Theorem 2.5 in Dawson and Li [S] are 
satished by taking the spaces; E = {1,2}, Uq = Ui = (0, oo] x (0, oo), the measures; 
7r(d^) = (5i(d2;) + (52(d2), pLQ{dz,du) = l{z<i}pL{dz)du, pLi{dz,du) = l{i<^<oo}/i(dx)dM + 
qSoo{dz)du, and the functions 

b{x) = a{x A n), a{x, z) =l{ 2 =i}'y 272 (x^\ri) + crl{ 2 = 2 }(a; A n), 

go{x, Z, u) = {ZA n)l{u<xAn}, 9l{x, Z, U) = {z A n)l{u<xAn}- 

Then, there exists a unique non-negative strong solution to ([9]) and therefore there exists 
a unique non-negative strong solution to ([8]). For m > 1 let Tm = inf{t > 0 : > m}. 

Since 0 < < m for 0 < f < r^, the trajectory t i-A zf”^^ has no jumps larger than 
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m on the interval [0,rm). Then, we have that satishes (j4]) for 0 < t < Tm- For 
n > m > 1, let (Yt, t > 0) be the strong solution to 


Yt =ZI^I + a (n A n)ds + / A 

Jo Jo 

a MYsAn)— ^ 

,3,1) Jo 


+ / (y; A n)dB2+, + 


+ 



2 ; A n)N{Tm + ds, dz, dn) 

(YsAn) — 

{z A n)N{Tm + ds, dz, dn). 


0 ^[l,oo]J0 


Next, we dehne for 0 < t < and for t > Tm- It is not difficult 

to see that is a solution to ([8]). By the strong uniqueness, we get that 

for all f > 0. In particular, Z^"^ = zf™^ < m, for 0 < f < Tm- Consequently, the sequence 
is non-decreasing. We dehne the process {Zt,t > 0) as follows 

* I oo ii t > lim Tm- 

m^oo 

Therefore, the process {Zt,t > 0) is a weak solution to (jl]). 

Now, we consider Z' and Z", two solutions of (j4]) and dehne 


Tm = inf{t > 0 : Zj' > m}, r" = inf{f > 0 : Z" > m}, 


and am = Tm/\Tm- Thus, Z' and Z" satisfy ([8]) on [0, dm), implying that both processes are 
indistinguishable on [0,crm)- If O’oo = hm am < oo, we have two situations to consider, 

m—)-oo 

either Z' or Z" explodes at (Too continuously or by a jump of inhnite size. If Z' or Z" 
explodes continuously then both processes explode at the same time since they are indis¬ 
tinguishable on [0, am), for all m > 1. If Z' or Z" explodes by a jump of inhnite size, then 
this jump comes from an atom of the Poisson random measure N, so that both processes 
have it. Since after this time both processes are equal to oo, and since the integral with 
respect to the Poisson process diverges, we get that Z' and Z" are indistinguishable. Then 
according to Theorem 137 in Situ [ 21 ], there is a unique strong solution to (jl]) that we 
denote by {Zt, t > 0). The strong Markov property follows since there is a strong solution, 
the integrators are Levy processes and the integrand functions are not time dependent (see 
for instance. Theorem V.32 in Protter [22|, where the Lipschitz property is just needed to 
guarantee the existence and uniqueness of the strong solution) and by Ito’s formula it is 
easy to show that the inhnitesimal generator of {Zt,t > 0) is given by fISl) . 

The branching property of Zj conditioned on K, is inherited from the branching pro¬ 
perty of the underlying CB-process. Its proof follows similar arguments as those used in 
Caballero et ah [7] (see p.77-79). More precisely, we consider (Z^, f > 0) and (Zj, f > 0) two 
solutions of dlj) starting from x and y, respectively, which are independent conditioned on 
the Brownian random environment > 0). In other words, we consider (5^)^ 

and the Brownian motions and Poisson random measures associated to Z and 

Z, respectively, which are mutually independent. Then conditioned on the environment. 



we observe 


:= 


Zt + Zt = x + y + a / Z',ds + / ) 

Jo Jo 

+ [ +a [ Z'.dBi^^ + Ut + Vt, 


where 


Ut : = 


r-^s- 


^A^*'^^(ds, dz, du) + 


Z.- 


ziV*'^^(ds, dz, dti) 


Jo J[l,oo]Jo Jo ~'[l,oo] Jo 

and Vt represents the sum of the integral terms with the compensated Poisson random 
measures. By the Levy’s Characterization Theorem, we deduce 


Wt = 


L{Z'^ 0 }- 


VW^s 


( 2 ) 


+ 


l{z'=o}d5l^) 


is a Brownian motion, since (hP)t = t, for t > 0. 

Next, we introduce M(ds,d 2 ,dM) an independent Poisson random measure, with in¬ 
tensity dsA(d^)dM and dehne 


7V(ds, d^;, dn) = l{„<z,_}iVd)(ds, d^;, dn) 

+ l{z,_<«<z'-}iV(^^(ds,d^,dM) -h l{z'_<^M(ds,d2;,dM). 


Using the same arguments as in [7] (p. 78-79), one can deduce that is a Poisson random 
measure with intensity dsA(d 2 ;)dM, and 

Ut = [ [ [ zN{ds, dz, dn) fo = 

Jo J[l,oo] Jo 

Therefore, given the environment, Z' satishes the SDE (jlj) with Zq = x + y, implying that 
it satishes the branching property. 

In order to deduce ([6]), we follow similar arguments as used in Bansaye [1]. To this 
purpose, we introduce Zf = Zte~^^ and take F G M+). An application of Ito’s 

formula guarantees that F{t, Zt) conditioned on iP is a focal martingale if and only if for 
every f > 0, 

0 = l‘ z.) + “Z. Af(<,, Z.) + Z,)) ds 

+ 11’°^- 2. + - F(s, Z.) - Af(s, A(d2)ds. 

If in addition, F is bounded, it will be a true martingale if the previous equality holds. 
By choosing F{s, x) = exp {—xvt{s, A, K)}, where Vt{s, A, K) is differentiable with respect 



fZ' 


zN{ds, dz, dn). 


0 J(o,i) Jo 
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to the variable s, non-negative and such that Vt(t, A, K) = A for all A > 0, we observe that 
(F(s, Zg), s < t) conditioned on is a martingale if and only if 


A, K) = - avt{s, A, K) + 7 ^(ni(s, A, K)Ye 

which is equivalent to Vt{s, A, K) satisfying ([7]). Since ^|J is locally Lipschitz on (0, oo), the 
existence and uniqueness of Vt follows from the Picard-Lindelof theorem. 

The above implies that the process {exp{—ZsVt{s, X, K)},0 < s < t) conditioned on K 
is a martingale, and hence 


E, 




E, 


exp 


-Zont(0, A, JP)| 


K 


exp{-znt(0, A, JP)}, 


which completes the proof. 


□ 


Remark 1. Observe that in the case when |'0'(O-|-)| < oo, the auxiliary process can be 
taken as follows 

+ mt, for t >0, 

where 

2 

m = -^/^'(0+) - y. 

Following the same arguments as in the last part of the proof of the previous Theorem 
and replacing K with one can deduce that Vt{s, A,iP®) is the unique solution to the 
backward differential equation 

yni(s, A, JP®) = e^"°VoK(s, A,iP(°))e"^™), (10) 

where 


fo{X) =-q +'y'^X^ + / (e - 1Aa;)/i(da;). 


(0,cxd) 


In this case, the process Z conditioned on satisfies, for every t > 0 and z > 0, 


E, 


exp 


XZtC 


-K 


( 0 ) 


K(o) 


= exp 


zvt(0,X,K(°^)j 


a.s. 


( 11 ) 


( 12 ) 


Before we continue with the exposition of this manuscript, we would like to provide 
some examples where we can compute explicitly the Laplace exponent of the CB-process 
in a Brownian random environment. 


Example 1. (Neveu case): The Neveu branching process in a Brownian random 
environment has branching mechanism given by ^{u) = ulog{u), for n > 0. In this 
particular case the backward differential equation ([7]) satishes 

yni(s. A, (5) = nt(s. A, 6) log{e~^’’vt{s, A, 5)). 
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Providing that Vt{t,X,6) = A, one can solve the above equation and after some straight¬ 
forward computation we deduce 


vt{s, A, (5) = exp <1 e"* ( / e “(5„dn log(A)e * 
Hence, from identity (jb]) we get 


E, 


exp j - AZte"^‘ 


K 


= exp < —zX^ exp s e ^Kgds 


a.s. 


(13) 


Observe that the r.v. e ^Kgds, is normal distributed with mean —^(1 — e ^ —te *) and 
variance ^(1 + 4e“* — 3e“^*), for t > 0. In other words, the Laplace transform of Zte~^^ 
can be determined by the Laplace transform of a log-normal distribution which we know 
exists but there is not an explicit form of it. 


Example 2. (Feller case): Assume that /i(0, cxo) = 0, thus the CB-process in a 
Brownian random environment dlj) is reduced to the following SDE 


Zt = Zo +a [ ZAs + a [ ZAB^^^ + [ 
Jo Jo Jo 

This SDE is equivalent to the strong solution of the SDE 


dZt =Y^tdt + ZAKt + ^/2AZABs, 
dKt =aodt -|- adBl:^\ 


where oq = a — cr^/2, which is the branching diffusion in random environment studied by 
Boinghoff and Hutzenthaler [6]. 

The backward differential equation ([^ satishes 

A, 6) = -avt{s, A, <5) -1- (s. A, 6)e~^\ 


If Vt{t, A, 6) = A, the above equation can be solved and after some computations one can 
deduce 


vtis, A, 6 ) = e-“" (^(Ae“*)-^ 


-1 


Hence, from identity (|6]) we get 


E. 


exp < — XZte 


— (Kt+at) 


K 


= exp 


ft \-i' 

f A"^ + 7 ^ / ) j> a.s. (14) 


The r.v. e is known as the exponential functional a Brownian motion with 

drift, and it has been deeply studied by many authors, see for instance HEmn]. 
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Example 3. (Stable case): For our last example, we assume that the branching 
mechanism is of the form 

'0(A) = —aX + , A > 0, 

for some {3 G (—1, 0) U (0,1), a G M, and is such that 

C;3<0 if0G(-l,O), 

c^>0 if0G(O,l). 

Under this assumption, the process Z satishes the following stochastic differential equation 


Zt =Zo + a I Z,ds + a / + 


t /*O0 pZs — 



zN{ds, dz, du) 


(15) 


0 Jo Jo 


where is a standard Brownian motion and 

jv(ds dz du) = I if /3 e (-1,0), 

’ ’ \ N{ds,dz,du) if/I e (0,1), 

where N is an independent Poisson random measure with intensity 

cii(3{(3 + 1 ) 1 


r(l-/ 3 ) z 2+/3 


dsdzdu. 


and N is its compensated version. 

In this case, we note 

0'(O+) = 


-oo if /3 G (—1, 0), 
-a if /3 G (0,1). 


Hence, when (3 G (0,1), we have = Kt + at, for f > 0. In both cases, we use the 
backward differential equation ([7]) and observe that it satishes 

A,(5) = -avt{s,X,6) + C/ 3 V^^^{s, X,6)e~^^’’. 
os 

Similarly to the Feller case, assuming that Vt(t, A, 5) = A, we can solve the above equation 
and after some straightforward computations, we get 

- 1//3 

Vt{s, A, (5) = e-“" ( (Ae“*)-^ + f3cp I 


- 1 // 3 - 


Hence, from ([6]) we get the following a.s. identity 


E, 


exp 


XZ^e 


-Kt 


K 


= exp < —z 


(^(Ae"*)-^ + (3cp 


which clearly implies the following a.s. identity 


E, 


exp < — XZtC 


— {Kt+at) 


K 


= exp 


(^A-^ + I3cfi 


- 1 // 3 - 


(16) 


We hnish this example by observing that the r.v. J^e 3 iKu+au)^^ jg exponential 
functional of the Brownian motion with drift {f3{Ku + au),u > 0). 
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3 Long-term behaviour 

Similarly to the CB-processes case, there are three events which are of immediate concern 
for the process Z, explosion, absorption and extinction. Recall that the event of explosion 
at hxed time t, is given by {Zt = cxo}. When < oo) = 1, for alH > 0 and 2; > 0, we 

say the process is conservative. In the second event, we observe from the dehnition of Z 
that if Zt = 0 for some f > 0, then Zt+s = 0 for all s > 0, which makes 0 an absorbing 
state. As Zt is to be thonght of as the size of a given population at time t, the event 
{limt^oo Zt = 0} is referred as extinction. 

Up to our knowledge, explosion has never been studied before for branching processes 
in random environment even in the discrete setting. Most of the results that appear in the 
literature are related to extinction. In this section, we hrst provide a sufficient condition 
under which the process Z is conservative and an example where we can determine ex¬ 
plicitly the probability of explosion. Under the condition that the process is conservative, 
we study the probability of extinction under the influence of the random environment. 

In our particular case, the events of explosion and absorption are not so easy to deduce 
in full generality. In the next section, we provide an example under which both events can 
be computed explicitly, as well as their asymptotic behaviour when time increases. 


3.1 Explosion and conservative processes 

Recall that V'XO-I-) G [—cxo, cxo), and that whenever |'0'(O-f-)| < 00 , we write 


m = —-^'(0-1-) 


a 


The following proposition provides necessary conditions under which the process Z is 
conservative. 

Proposition 1. Assume that g = 0 and |'0XO-f-)| < 00 , then a CBBRE with branching 
mechanism 'ip is conservative. 

Proof. Recall that under our assumption, the auxiliary process takes the form 


+ mf, for f > 0, 


and that Vt{s, X, is the unique solution to the backward differential equation flTOj) 
with q = 0. From identity (IT^ . we know that 




exp 


XZtC 


-K, 


(0) 


= E 


exp 


zvt{^, A, 


Thus if we take limits as A 1 0, we deduce that for 2; > 0 


Fz(yZt < oo'j = limE^ 
A4..O 


exp < — XZtC 


-K, 


(0) 


= E 


exp u(0 






where the limits are justihed by Monotonicity and Dominated Convergence. This implies 
that a CBBRE is conservative if and only if 

limnt(0,A,A:®) = 0. 

A4.0 
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Let us introduce the function $(A) = A“^'0o(A), where -00 is given by ffTT]) and observe 
that <h(0) = -00(0+) = 0- Since 'ipo is convex, we deduce that <h is increasing. Finally, if 
we solve equation (IT0|) with '0o(A) = A<F(A), we get 

nt(s, A,iF) = A exp |- <F(e"^"nt(r, A,iF))dr|. 

Therefore, since <h is increasing and *h(0) = 0, we have 

0 < limnt(0. A, JF) = limAexp/— [ <F(e“^’'nt(r, A, JF))drl < limA = 0, 

A—^0 A—>-0 ^ ^ Q j A—>-0 

implying that Z is conservative. 

□ 


Recall that in the case when there is no random environment, i.e. a = 0, we know that 
a CB-process with branching mechanism ^|J is conservative if and only if 

f du _ 

In the case when the random environment is present, it is not so clear how to get a 
necessary and sufficient condition in terms of the branching mechanism since the random 
environment is affecting the monotonicity of ijj in the backward differential equation ([7]). 

We now provide two interesting examples in the case when '0'(O+) = —oo, that behave 
completely differently. 


1. Stable case with jS G (—1,0). Recall that in this case = —au + 
where a G M and cp is a negative constant. From straightforward computations, we get 


'0'(O+) = -cx). 


and 


du 


< oo, 


Jo+ \Hu)\ 

On the other hand, from identity flTBD and taking limits as A J, 0, we deduce that for z > 0 

rt \ - 1 // 3 ' 


(^Zt < OO 


K ) = exp < —z 


implying 


FJ Zt = oo 


K ) = 1 — exp 




a.s., 


(17) 


- 1 // 3 - 


> 0 . 


In other words the stable CBBRE with fd G (—1,0) explodes with positive probability 
for any t > 0. Moreover, if the process (iF„ + (y.u,u > 0) does not drift to —oo, i.e. 
a — cr^/2 > 0, we deduce from Theorem 1 in Bertoin and Yor [1], 

g-/3(x„+Q«)^^ —)■ CX), as t ^ oo. 
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implying 


a.s. 


lim Zt = CX3, 

t—^OO 

On the other hand, if the process {Ku + au,u > 0) drifts to —oo, i.e. a — < 0, we 

have an interesting long-term behaviour of the process Z. In fact, we deduce from the 
Dominated Convergence Theorem 


P 


2: 




1 -E 


exp 




r*00 \ ~1//^ 


The above probability is positive since 

g-/3(-K'„+aw)^^ < OO a.s., 

according to Theorem 1 in Bertoin and Yor |1]. In this particular case, we will discuss the 
asymptotic behaviour of the probability of explosion in Section 4. 



2. Neveu case. In this case, recall that 'ip{u) = ulogu. In particular 


'0'(O-f-) = —00 and 



du 

u logu 


= 00. 


By taking limits as A 0 in ffT3|) . one can see that the process is conservative conditionally 
on the environment, i.e. 


P^ {Zt < oo|iC) = 1, 


for all t G (0, 00) and z G [0, 00). 


3.2 Extinction probabilities 

Here, we consider CBBREs that are conservative. The following result provides a criteria 
that depends on the behaviour of the auxiliary process which allows us to compute 
the probability of extinction of a CBBRE. Recall that the event of extinction is equal to 
{hmt^^Zt = 0 }. 


Proposition 2. Assume that g = 0 and |'0'(O-|-)| < 00 . Let {Zt,t >0) be a CBBRE with 
branching mechanism given by ijj and z > 0. 

i) // m < 0, then P^ { lim = 0 = 1, a.s. 

\t^OO J 

ii) //m = 0, then P^ (liminfZt = 0 ] = 1, a.s. Moreover if 'j > 0 then 

\ t^OO J 


P^ lim Zi 

\t—^oo 



1, a.s. 
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in) // m > 0 and 


a:ln(a:) ^{dx) < oo, 


then liminfZi > 0 

\ t^OO 

such that 


i^(o) 

> 0 a.s., and there exists a non-negative finite r.v. W 


Ztc 


-K, 


(0) 


t—>00 


a.s and {W = O} = | Um Zt = o|. 


Moreover if j > 0, we have 


P2 (lim Zt = 0 > 0, a.s. 

\t^OO J 


and, 


lim Zt = 0 ) > I 1 + 

\t^OO 


za 


7 


2 \ -- 


Proof. Recall that under our assumption, i.e. |'0XO+)| < 00 , the auxiliary process can be 
written as 

+ mt, for t > 0, 

and the function Vt{s, A, satishes the backward differential equation flTOl) . 

Similarly to the last part of the proof of Theorem 1, one can prove that ZtC * is 
a non-negative local martingale. Therefore Zte~^*- is a non-negative supermartingale 
and it converges a.s. to a non-negative hnite random variable, here denoted by W. This 
implies the statement of part (i) and the hrst statement of part (ii). 

In order to prove the second statement of part (ii), we observe that if 7 > 0, then the 
backward differential equation fllOl) satishes 


A,^®) > 7^nt(s, A,R'®)^e \ 

os 


Therefore 


Vt 


which implies the following inequality, 

FfiZt = 0|R:(°)) > exp I^ e-^^°^ds^ 
Since m < 0, we deduce from Theorem 1 in Bertoin and Yor [1] 


-1 


e ds —)■ cxo, as f —)■ cxo. 


(0) 


(18) 


( 19 ) 


In consequence, we have 


P, (lim Zt = 0 RT®) = 1 

\t—>-oo / 


a.s. 
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Now, we prove part (iii). We first note that Vt{-, A, the solution to the backward 

differential equation (1T0|) . is non-decreasing on [0,t] (since 'ipo is positive). Thus for all 
s E [0,t], nt(s, < A. From the proof of part (ii) in Proposition 1, we know that 

the function ‘F(A) = A“^'0o(A), is increasing. Hence 

A, = vt{s, A, A, < Vt{s, A, JF(°))<h(Ae-'^«“'). 

os 

Therefore, for every s <t, we have 


Vt{s, X, > A exp 



In particular. 


lim inf (0, A, > A 


t^OO 


exp 


<F(Ae ^)ds 


According to Proposition 3.3 in Salminen and Yor 1211 , we have 


<h(Ae )ds < cxo a.s., if and only if 


(0) 


$(Ae ^)dy < oo. 


Observe 


<h(Ae ^)dy = 


m 


de 


= 7'A + 


de 


0 ^ J (0,oo) 


(e — 1 + 9x)y.{dx) 


= 7^A -I- / fi{dx) I {e — 1 + 9x) 


= 7'A + 


(0,oo) 

f X 

(0,oo) \J0 


¥ 


(e ^-1 + y)^ j/u(dx). 




Since the function 


9x{x) = [ (e ^ - l + y)^, 

Jo r 


is equivalent to Aa;/2 as x —)■ 0 and equivalent to Inx as a: ^ cxo, we deduce that 


/ <h(Ae ^)dy < oo if and only if 

Jo 

In other words, 

/ <h(Ae )ds < oo a.s., if and only if 

Jo 

If the integral condition from above is satished, then 
hminfnt(0,A,A:(°)) > A exp 


X In x/i(dx) < oo. 


X In xfi{dx) < oo. 


t—^OO 


$(Ae-^“')ds > 0, 


(0) 
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implying 


E, 


-,-xw 


kw 


< exp < —z A exp 


<l>(Ae-^«“')ds < 1 , 


( 0 ) 


and in particular P^lliminfZi > OiPW > 0 a.s. Next, we use Lemma 20 in [1] and the 
V t^oo / 

branching property of Z, to deduce 

{PL = 0} = = o|. 

Now assume 7 > 0. From inequality (1T8|1 and the Dominate Convergence Theorem, we 
deduce 


(lim Zt = 0 > 

\t—^oo J 


> exp < —z I 7 




-1 


a.s. 


On the other hand, since m > 0, we deduce from Theorem 1 in Bertoin and Yor |1] (or 
Proposition 3.3 in Salminen and Yor 


-.-K, 


(0) 


ds < cxo, a.s.. 


( 20 ) 


implying that P2(hmi^oo Zt = > 0, a.s. Finally, according to Dufresne [9], 


.-K 


(0) 


ds has the same law as (2F 2 


-1 


( 21 ) 


where F.^ is Gamma r.v. with parameter v. After straightforward computations, we deduce 


P, ( lim Zi- = 0] > ( 1 + 


\t^OO 


Z(J^ 


7 


The proof of Proposition 2 is now complete. 


□ 


Remark 2. When m > 0, we also have an upper bound for the probability of extinction 
under the assumption 


7 > 0 and K = x‘^p,{dx) < 00. 


More precisely, under the above assumption we have, for z > 0 


1 + 


za 


27^ J Vt—>00 

The upper bound follows from the a.s. inequality 

d 


< P2 (" lim Zi = 0) < (1 + 


2 7^ + K 


ds 


Vt{s, A, < (7^ + K)vt{s, A, 


-K 


(0) 


2m 
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It is important to note that when the branching mechanism is of the form 'i/j{u) = 
—au + for {3 G (0,1], one can deduce directly from ffl^ and flT^ by taking A and 

f to oo, and (fT9D that 

lim Zt = 0, a.s., for m < 0. 

)-oo 

Similarly, using fl20|) . in the case when m > 0 we have 

( ^oo \ —1//3 'i 

J L a.s., 

and in particular 

P = 0j = P^ ^ lim Zt = 0^ = 

The latter probability can be computed explicitly using ([21]). 

We hnish this section with a remark on the Neveu case. If we take limits as A t cxo, 
in flTS]) . we obtain that the Neveu CBBRE survives conditionally on the environment, in 
other words 



P2 ( lim Zt = Ol/’i') = 
\t^oo J 


exp 


(^t > o|-^) — 1) 

for all t G (0, cxo) and 2; G (0, cxo). Moreover since the process has cadlag paths, we deduce 
the Neveu CBBRE survives a.s., i.e. 

^zi^Zt > 0, for all t > oj = 1, for all z > 0. 

On the one hand, using integration by parts we obtain 


e-^K,ds = a I - (1 - - e-^Kt 

7n 7 


Since 


a 


e ‘dBt’ = y (1 


we have that 


lim 

t^OO 


e ^0 < 00 


e ^KAs = a 


and 


E 


a 


^ 2 


ABi^^ 


cr" 

T’ 


2 2 

exists and its law is Gaussian with mean — ^ and variance Hence, if we take limits as 


f t 00 (II3|), we observe 


E. 


exp s — A lim Ztc 


-K, 


t^OO 


K 


= exp < —2; exp 


e ^KAs 


for all z > 0. 


Since the right-hand side of the above identity does not depend on A, this implies that 


P2 ( lim Zte ^* = 0 

. t—>00 


K) = exp < —2; exp 


e ^iC,ds 


for all z > 0, 
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for all z > 0. 


and taking expectations in the above identity, we deduce 


P. 




E 


exp < —z exp 


^JCds 


In conclusion, the Neveu process is conservative and survives a.s., but the extinction 
probability is given by the Laplace transform of a log-normal distribution. 

Finally, if we multiply flT^ by e^‘, differentiate with respect to A and then we take 
expectations on both sides, we deduce 


E, 


Zt exp I - XZte | 
= ze~^X 


^E 


exp ^ Kt + J e ‘^Kgds — z\^ exp 


e-^KAs 


Taking limits as A goes to 0, it is clear 

E2 [Zt] = 00 , for t > 0, z > 0. 


4 Stable case 


The stable case is perhaps one of the most interesting examples of CB-processes. One of 
the advantages of this class of CB-processes is that we can perform explicit computations 
of many functionals, see for instance [aittniE], and that they appear in many other areas 
of probability such as coalescent theory, fragmentation theory. Levy trees and self-similar 
Markov process to name but a few. As we will see below, we can also perform a lot 
of explicit computations when the stable CB-process is affected by a Brownian random 
environment. 

In the sequel, we shall assume that the branching mechanism satishes 'ip{u) = —au + 
Ci 3 U^^^, for M > 0, and f3 E (—1, 0) U (0,1). Recall that a G M, and 


f C/3 < 0 if j3 E (-1,0), 

I C/3 > 0 if /9 e (0,1). 

Boinghoff and Hutzenthaler [6] studied the particular case when /3 = 1, also known as 
the Feller diffusion case. The authors in [6] gave a precise asymptotic behaviour for the 
survival probability and also studied the so-called Q-process. In this section, we prove 
similar results for the case when /3 E (0,1) and we obtain new results on the asymptotic 
behaviour of non-explosion for the case when [5 E (—1, 0). 

Recall from identity ffTOj) that the stable CBBRE Z = {Zt,t > 0), satisfies 


E, 


exp < — XZte 


— {Kt+at) 


K 


= exp 


( A"^ -f 


u 


'0 


-1//3- 


If we take limits as A goes to 00, in the above identity we obtain for all z, t > 0, 

rt 


F,{Zt>0\K) =1- exp I 


L{/3>0}, 


a.s., (22) 


20 











where the same holds true for the Feller case by taking (3 = 1 and C/3 = 7^. 

On the other hand, if we take limits as A goes to 0, we deduce that t, z > 0 


i^Zt < 00 


K ) = exp 




-1//3' 


l{/3<0} + l{/3>0} a.S. 


It is then clear that if /3 G (—1, 0), then the survival probability equals 1, for all f > 0. If 
(3 G (0,1] then the process is conservative. 

From the above identities, a natural question arises: can we determine the asymptotic 
behaviour, when t goes to 00 , of > 0) for {3 G (0,1] and < 00) for (3 G 

(—1,0)? We observe below that the answer of this question depends on a hne study of 
the asymptotic behaviour of 

\ - 1//3 

For this purpose, let us recall some interesting facts of the exponential functional of a 
Brownian motion with drift. 



E 


exp 


(3cp 


4.1 Exponential functional of a Brownian motion with drift 


In what follows, the following functional will be of particular interest. Let be the 
exponential functional of a Brownian motion with drift 7 G M, in other words 


j(v) _ 

-'t •— 



s + Bs) Us, 


t G [0, cxo). 


The law of such random variable have been deeply studied by many authors. Up to our 
knowledge this is the unique example for which there exist an explicit formula for the joint 
distribution of + rjt), see for instance Proposition 2 in Matsumoto and Yor [T9] . 

In particular, for all t, u E (0, 00) and x G M, we have 


P (E du\Bt + rjt = x] = 


y/2nt 


u 


X 

exp ^ exp 




2x 


0e^/u{t)du, 


where 


6r{t) = , ^ 6 2* [ 6 2* ^“®*^*^^Uinh(|/) sin ) dy, r > 0. 




\ t J 


The following lemma generalizes Lemma 4.2 in [6] from p = 1 to p > 0. 
Lemma 1. Let p G M and p > 0. Then for every t > 0, we have 




ii 


E 

E 


riv) 


riv) 


-p 


-2p 



l^j(-(n-2p))^ P 

< g(2p^-2p??)t£ 



E 


r((»7-2p)) 

n/2 


-P 
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Proof. Using the time reversal property for Brownian motion, we observe that the process 
{rjt + Bt — r}{t — s) — Bt-s, 0 < s <t) has the same law as {rjs + Bg, 0 < s < t). Then, we 
deduce that 


g 2 (? 7 s+Bs)^g has the same law as / e ds. 


We now introduce the exponential change of measure known as the Esscher transform or 
Girsanov’s formula 


dp(^) 

dP 



for A e M, 


where {Bt)t>o is the natural hltration generated by the Brownian motion B which is 
naturally completed. Observe that under pG)^ the process i? is a Brownian motion with 
drift A. Hence, taking A = —2p, we deduce 


E 


= E 

^-2p(rit+Bt) f f g-20s+Ss)^g^ 


\Jo J 


\Jo ) 


= g-2p»7ig2p2tjg(_2p) 






g-2((p-2p)s+Bs)^g 


-P 


-P 


which implies the hrst identity, thanks to the symmetry property of Brownian motion. 
In order to get the second identity, we observe 


rt/2 


ft/2 


Jo Jo Jo 

where Bg = Bgj^t /2 — -Bt/ 2 , s > 0, is a Brownian motion which is independent of (!?„, 0 < 
u < t/2). Therefore, using part (i), we deduce 


E 


-t \ - 2 p' 


'0 


< E 


rt/2 \ P' 

gPt+ 2 St /2 / Ms 


E 




ft/2 


gPt+ 2 Bi /2 / 


\il 

rt/2 


-P- 


0 

y 


-p- 





E 

^j(-0-2p))^ 

r 


On the other hand from the Esscher transform with A = —2p, we get 
E 


rt/2 \ P' 

gPt+ 2 St /2 / 


= e-PP^eP^^E^-^P^ 
= e-PP^eP^^E 


t/2 \ P 

g2(ps+S.)ds 


»t/2 


-PI 


g2((p-2p)s+Bs)^g 
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Putting all the pieces together, we deduce 


E 


ft \ -V 




j({r)-2p)) 

^tl2 


-P 


E 


t{-{v-2p)) 

-'t/2 


-P 


The proof of the lemma is now complete. 
We are also interested in 


□ 


/(»?) 


exp 


s + Bs) Us, 


which is hnite a.s. whenever rj < 0. We recall that according to Dufresne [9], 

has the same law as (2r_,j) , 


(23) 


where P^ is Gamma r.v. with parameter v. 


4.2 Explosion probability 

Throughout this section, we assume that [3 G (—1,0). As we will see in the main result 
of this section, the asymptotic behaviour of the probability of explosion depends on the 
value of 

m = a-. 

2 

We recall that when a = 0, i.e there is no random environment, the stable CB-processes 
explodes with positive probability. Moreover, when a = 0, we can compute explicitly the 
asymptotic behaviour of the probability of explosion. 

When a Brownian random environment affects the stable CB-process, the process be¬ 
haves completely different. In fact, it also explodes with positive probability but we have 
three different regimes of the asymptotic behaviour of the non-explosion probability that 
depends on the parameters of the random environment. Up to our knowledge, this be¬ 
haviour was never observed or studied before. We call these regimes suhcritical-explosion, 
critical-explosion or supercritical-explosion depending on whether this probability stays 
positive, converge to zero polynomially fast or converges to zero exponentially fast. 

Let 


r] := 


(3a‘^ 


-m 


and 


k = 


2c[- 


2 \ 1//9 


and dehne 


5f(a:) := exp , for x > 0. 

From identity flTT)) and the scaling property, we deduce for f3 G (—1,0) and ?7 > — 1, 


P2 {Zt < 00 ) = E 
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(v) 




g{z^v)ptaii^,p{v)^v, 


(24) 
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where denotes the density function of l/2llP'^ which according to Matsunioto and Yor 
[T9] , satishes 




g-r,2^/2g7rV2t^ / r] + 2 


V2' 






g?V2i'g(f?-l)/2g-3;s 


(25) 


'0 JO 


sinh(0cosh(0sin(7r^/z/)^^^^ 

^ / I , /^^2^2±2 ClsClS. 

(s + cosh(4)^) 2 


We also denote 




E 



for 6* > 0. 


Theorem 2. Let {Z+A > 0) he the stable CBBRE with index B G (—1,0) defined by the 
SDE ([ID with Zo = z>0. 


i) Subcritical-explosion. //m<0, then 


liniP^ {Zt < oo) 

t^OO 




a) Critical-explosion. //m = 0, then 


vTr/3(T Jo X ■ 

Hi) Supercritical-explosion. //m>0, then 

3 m^f 8 / R \ 

limf2ei^P^ {Zt <oo) = / g{z^v)(t)j^{,v)<lv, 

t->-oo p^a Jo 


lim y/t {Zt < oo) = 

t^OO 


where 

Brjiv) 




ooshg); 
{u + v cosh(,^)2)'’2 


Proof. Our arguments follows from similar reasoning as in the proof of Theorem 1.1 in 
Boinghoff and Hutzenthaler [6]. For this reason, following the same notation as in [6], we 
just provide the fundamental ideas of the proof. 

The subcritical-explosion case (i) follows from the identity in law by Dufresne fl2^ . 
More precisely, from fl23|) . fl2T|) and the Dominated Convergence Theorem, we deduce 


limP^ {Zt < oo) = E 


t^OO 


exp { -zkVjJ 


= Crt,0 {zk) 


In order to prove the critical-explosion case (ii), we use Lemma 4.4 in [6]. From identity 
and applying Lemma 4.4 in [6] to 


g{z^x) = exp {—< 


X 


-1//3 


Zk ’ 


for x > 0, 


(26) 


24 















we get 


lim Vt P2 (Zt < 00 ) 

>00 


2 

hm \ - 

f3a V 4 


E 




13 


orW 

^-'/32o-2t/4, 


y/nf3a Jo X ’ 


which is finite since the inequality fl26|) holds. 

We now consider the supercritical-explosion case (iii). Observe that for all n > 0, 


g{z^x) = exp {—for a: > 0. 


Therefore using the above inequality for a fixed n, Lemma 4.5 in [6] and identity fl24l) . we 
obtain that for each 0 < m that satisfies m < ?7,(T^/2, the following limit holds 


limt3/2 ^ 

t—>00 


limf3/2 

t^OO 


9 


I 21^^^ / 


8 





where (j)rj is defined as in the statement of the Theorem. Since this limit holds for any 
n > 1, we deduce that it must hold for m > 0. This completes the proof. 


□ 


4.3 Survival probability 

Throughout this section, we assume that (3 G (0,1). One of the aims of this section is 
to compute the asymptotic behaviour of the survival probability and we will see that it 
depends on the value of m. We find five different regimes as in the Feller case (see for 
instance Theorem 1.1 in [B]) and CB-processes with catastrophes (see for instance Propo¬ 
sition 5 in m)- Recall that in the classical theory of branching processes, the survival 
probability stays positive, converges to zero polynomially fast or converges to zero expo¬ 
nentially fast, depending of whether the process is supercritical (m > 0), critical (m = 0) 
or subcritical (m < 0), respectively. In the stable CBBRE there is another phase tran¬ 
sition in the subcritical regime. This phase transition occurs when m = —We say 
that the stable CBBRE is weakly subcritical if —< m < 0, intermediately subcritical if 
m = —cr^ and strongly subcritical if m < — 

Recall that 


and define 



and 


k = 


2cp) 


/(x) := 1 — exp {—kx^/^} , for x > 0. 
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From identity fl^ and the scaling property, we deduce for /? > 0 and 77 > —1, 


(^t > 0) 


E 


/ 


{ ] 

\ ] 

\‘^-'/32<72i/4/ J 




( 27 ) 


where py^r^ denotes the density function of 1/21^'^ and is given in fl 2 ^ . 

Theorem 3. Let {Zt,t > 0) be the stable CBBRE with index fi G (0,1) defined by the 
SDE / fiTl) with Zq = z > 0. 

i) Supercritical, //m > 0, then 


limP.(Z,>0) = l-f;h£!^h|-^, 

t^Qo n\ 1 (— 77 ) 


n=0 


(28) 


a) Critical. //m = 0, then 


lim \/t P 2 {Zt > 0 ) = 


\/2 >;^(-zk)^ 


t^OO 


:E 


0r/3(j^ n\ \fi 

^ n=l ^ 


n? 


(29) 


Hi) Weakly subcritical. J/mG then 


hmf2ef^P,, {Zt > 0) = —^ 

t^oo 


f{z^v)(pr,{v)dv, 


where 


(30) 


POD /*OD 


(V = 


'0 Jo 


V27r 


-r 


^i«h(0 c°sh(0? 

2 / (u + 7;cosh(^)2)’’2 


iv) Intermediately subcritical. Ifm=—a‘^, then 

lim (Zt > 0) = z 

t—^OO 


V2 


kF . 


(31) 


v) Strongly subcritical. If m < —then 


lime-5(2-+-^hp^ (Zt > 0) = zk ^— 

^ ^ r{p-2/fi) 


t—^oo 


(32) 


Proof. As in the proof of Theorem O and following the same notation as in [ 6 ], we just 
provide the fundamental ideas of the proof. 

The supercritical case (i) follows from the identity in law by Dufresne (1231) . More 
precisely, from fl23|) . fl27)) and the Dominated Convergence Theorem, we deduce 


hmP^ (Zt > 0) = E 1 — exp zkF_,^| = 1 — y~~^(—^k) 


-yhV^/h\ = 


n=0 


. ni-v) 

77 .!F(— 77 ) ' 
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In order to prove the critical case (ii), we use Lemma 4.4 in [B]. From identity flTH) and 
applying Lemma 4.4 in [6] to 

g{x) := 1 — exp {—, x > 0, (33) 


we get 


lim \/i P 2 {Zt > 0) 

t^OO 


2 ii,„ E 

jSa V 4 


riv) 


V2 

\/^/Scr 


exp 


1 - exp |-zk 

{ j ^—dx. 


-1//3 


By Fubini’s theorem, it is easy to show that, for all g > 0 




n=l 





which implies (129|) . 

We now consider the weakly subcritical case (hi). Recall that inequality (l33|) still holds, 
then using Lemma 4.5 in [B] and identity (127)) . we obtain 


limt^/^ {Zt > 0) 

)-oo 


hmt3/2 

t^OO 


9 


I 2&^ / 

\^-'/32cr2t/4/ J 


o noo 

Jo ~ <P^{v)dv, 


where (prj is dehned as in the statement of the Theorem. 

In the remaining two cases we will use Lemma 4.1 in [6] and Lemma [H For the 
intermediately subcritical case (iv), we observe that rj = 2/f3. Hence, applying Lemma [T] 
with p = 1//9, we get 


E 




and 


E 



-2//3' 


- TTrt-n 


< e ^ E 



Now, applying Lemma 4.4 in [B], we deduce 


lim Vi E 

t—)-oo 




1 

_ _ = 

V271 a 




and 


lim Vi E 

t^OO 




0 . 


Therefore, we can apply Lemma 4.1 in [6] with q 


Vi and F) 



and 
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obtain 


(Z* > 0) = 


t—^OO 


£—)-oo 


1 - exp -zk (21 


r(v) \ 


= lim 

£—)-oo 


- ^ vr 

V/5 


21 


iv) 




-1//3- 


Finally for the strongly snbcritical case, we nse again Lemma [U with p = l/f3. First 
observe that p — 2/(3 > 0. Thns, the Monotone Convergence Theorem and the identity of 
Dnfresne ([23]) yield 


limE 

t —>00 


2L 


{-(7?-2//3)) 


t/2 


-IIP 


= E 


(2/( 


-(r,-2//3))j-V/3' 


E[(r,-2/^)‘"] 


r(i; - i/H) 
r(r,-2//J)' 


Since Ip goes to oo as t increases, from the Monotone Convergence Theorem, we get 

- 1//31 


limE 

t—IOO 




= 0 . 


Hence by applying Lemma 4.1 in [6] with ct = e Yt = (2J, 


obtain that 


riv) 


-i/p 


we 


t^OO 


lime 2(2m+(T2)tp > Q) = limCfl2o.2i/4E 

' t^oo ^ ' 

- 1//3) 

nv- 2 /A 


1 - exp -k 21 


riv) 


',320-24/4 


-1/p 


This completes the proof. □ 

4.4 Conditioned stable CBBRE 

Here, we are interested in stndying two conditioned versions of the stable CBBRE; the 
process conditioned to be never extinct (or Q-process) and the process conditioned on 
eventnal extinction. Onr methodology follows similar argnments as those nsed in Lambert 
ini and extend the resnlts obtained by Boinghoff and Hntzenthaler [6] and Hntzenthaler 
[12] in the stable case with f3 G (0,1). In particular, we obtain that the supercritical 
stable CBBRE conditioned on eventual extinction possesses a similar phase transition as 
the snbcritical stable CBBRE. It is important to note that such a phase transition has 
not been reported in the discrete case. For the continuous case, it was only observed in 
[T2] for (3 = 1. In contrast with the snbcritical regime, the phase transition is given at 
We say that the supercritical stable CBBRE conditioned on eventual extinction is 
weakly supercritical if 0 < m < f3a'^, intermediately supercritical if m = f3a'^, and strongly 
supercritical if m > /dcr^. 

















4.4.1 The process conditioned to be never extinct 

In order to study the stable CBBRE conditioned to be never extinct, we need the following 
Lemma. 


Lemma 2. For every t > 0, Zt is integrable. 

Proof. Differentiating the Laplace transform flT^ of Zt in A and taking limits as A goes to 
0, on both sides, we deduce 

E, [Zt\K] = 

which is an integrable random variable. □ 


Recall that 


T] = 


(3a‘ 


-m and k = 


/dcr^ 

2cfi 


2 \ 1//3 


We now dehne the function U : [0, oo) (0, oo) as follows 


V 2 f n 


0F/3( 




n=l 

f*oo 


n\ 




if m = 0, 


U(z) = 




kr 


zk 


V2 

0r/3(T 

r(r/ - 1//3) 


nv - 2 // 9 ) 


(j)n{v)dv if m G (—cr^, 0), 
if m = —cr^, 
if m < —cr^. 


where the function (/>„ is given as in Theorem [3l We also introduce 


0 


m 


6 := 6'(m, cr) = ■^ 

2m + 


if m = 0, 
if m G (-cr^O), 

if m < —a^. 


Let {iFt)t>o be the natural filtration generated by Z and Tq = inf{t > 0 : Zj = 0} be 
the extinction time of the process Z. The next proposition states, in the critical and 
subcritical cases, the existence of the Q-process. 


Proposition 3. Let {Zt,t > 0) be the stable CBBRE with index f3 G (0,1) defined by the 
SDE [151) with Zq = z > 0. Then for m < 0.- 

i) The conditional laws (• | Tq > t + s) converge as s ^ oo to a limit denoted by P^, 
in the sense that for any t >0 and A G iFt, 


lim P^ (A I To > t + s) = P^ (A). 

s—)-oo 
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a) 


The probability measure can be expressed as an h-transform of P based on the 
martingale 

Dt = e^^U{Zt), 


in the sense that 



Dt 

Uiz) 


dP. 


X't' 


Proof. We first prove part (i). Let z,s,t > 0, and A G Dt- From the Markov property, we 
observe 


Pz (A I To > t + s) 


P^ ^A n {To > t + 
Pz (^0 > t + s) 




Fzt {Zs > 0 ) ^ 

P. (Z,+. > 0) * ■ 


(34) 


On the other hand, since the mapping t i-G> is increasing and the function f{x) = 
1 — exp I —\<ix^D j ig decreasing, we deduce from IHT} and the Markov property that for 
any z,y > 0, 


0 < 


Py {Zs > 0) 

i^t+s > 0) 


E 


1 - exp -yk 21 


riv) 


2/32s/4 


- 1//3 


E 


1 - exp < -zk (21 


riv) 


- 1//3 


< 


ykE 

f2&"’ ^ 


E 

1 — exp < 

-zk (^ 27 ^ 2 ) 32 ^/ 4 ^ 

-V/i'l 1 

IJ 


Moreover, since diverge as t goes to oo, we have 
-i/M 


zkE 




2I_2o2 


< E 


1 - exp { -zk (21 


riv) 


2/32s/4 , 





< zkE 


21 


iv) 




-1/A 


for s sufficiently large. Then for any s greater than some bound chosen independently of 
Zi(a;), we necessarily have 

Pz.(Z.>0) ^2 
- P. (Z,+. > 0) - 2 " 


Now, from the asymptotic behaviour 

lim 


, (l3B and (l32D . we get 


Fz, (Z, > 0 ) e^*U(Zt) 


s^ooP^ (Zt+s > 0) [/(z) 

Hence, Dominated Convergence and identity flM|) imply 

e^^U(Zt), 


lim P 2 (A I To > t + s) = E^ 


L D(z) 


(35) 
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Next, we prove part (ii). In order to do so, we use fl5^ with A = to deduce 

E, [e^^U{Zt)] = U{z). 


Therefore, from the Markov property, we obtain 


E, 


J-J =e^^Ez. [e^^U{Zt)] =e^^U{Z,), 


implying that D is a martingale. 


□ 


4.4.2 The process conditioned on eventnal extinction 

Here, we assume that m > 0 and dehne for z > 0. 

TT , \ _ ^{-z]i.YT{n/(5 -T]) 

ri\ T{-r]) ■ 

n=0 ^ 

In the supercritical case, we are interested in the process conditioned on eventual extinc¬ 
tion. 

Proposition 4. Let {Zt,t > 0) be the stable CBBRE with index f3 G (0,1) defined by the 
SDE [TR] with Zq = z > 0. Then for m > 0, the conditional law 


P:(-) = P, (• I To < oo) , 


satisfies for any t >0, 


dPtl = ^li^dP. 


Ufiz) 

Moreover, {UfiZt),t >t)) is a martingale. 

Proof. Let z, t >0 and A G T*, then 


PjA 


P, (A n {To < CX)}) P, (A n {Zi+, = 0}) 

Jo < oo =-;-^-= hm---, 

' P,(To<oo) ™ Ufiz) 


On the other hand, the Markov property implies 


Pz (A n {Zt+s — 0}) — E 2 


Pz ( Zt+s — 0 


*Tt) 1 a 


= E. 


PzT^ = 0 )l 


Therefore using the Dominated Convergence Theorem, we deduce 


Pz A 


E, 


To < 00 I = hm 


FzAZs = 0]1 


E, 


U*{Zt)lA 


Ufiz) 


Ufiz) 


The proof that {U^{Zt),t > 0) is a martingale follows from the same argument as in the 
proof of part (ii) of Proposition |3l □ 
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Observe that F*(Zt > 0) goes to 0 as t —)■ cxo. Hence a natural problem to study is the 
rates of convergence of the survival probability of the CBBRE conditioned on eventual 
extinction. As we mentioned before, we obtain a phase transition which is similar to the 
subcritical regime. 

It is important to note that the arguments that we will use below also provides the 
rate of convergence of the inverse of exponential functionals of a Brownian motion with 
drift towards its limit, the Gamma random variable. The latter comes from the following 
observation. Since U^,{Zt) is a martingale, we deduce 


Fl{Zt > 0) =P^ (Zt > 0 I To < cx)) = 
1 


U,{z) 


E 


exp <! } 


U,{z) 
-E 


{U,{z)-F,{Z, = 0)) 


exp <; -zk ( 


- 1//3 


Another important identity that we will use in our arguments is the following identity in 
law 


r(^) 


,*>o S 


r(-'n) 


2r_^,f > 0 


where r_^ and ij: are independent (see for instance identity (1.1) in Matsumoto and 

Yor [I9]). We also introduce 


h{x, y) = exp {—— exp |—k {x + , 


x,y>Q. 


Then 


K{Zt > 0 ) 



h 





0 


(-»?) 


2 T 2 o 


(36) 


where r_^ and l\ are independent. 

Theorem 4. Let {Zt,t > 0) be the supercritieal stable CBBRE with index (3 G (0,1) 
defined by the SDE / fI3]) with Zq = z > 0. 

i) Weakly supercritical. If m E then 


limf2e^P*(Zi > 0) 

t—)-oo 


8 

/3^aW{\y\)Ufiz) 




h{z^x,z^y)(j)\r,\{y)x ^dxd?/. 


where 




sinh(0 cosh(0e 

(u + |/cosh(^)2)“ 


d.^dn. 


a) Intermediately supercritical. Ifm = (3a‘^, then 


t—^OO 


zk^/2 -^ (-^k)^ / n + 1 \ 

n\ \ fi J 
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in) Strongly supercritical. Ifm> (3a^, then 


t^OO 


-zk(rj + 2} ^ (-zk)” 
/3C/.(2)r(->))^ n! 


fn + l 


-p- 


1 


Proof. Similarly as in the proof of Theorems [2] and [3l and following the same notation as 
in [6], we just provide the fundamental ideas of the proof. 

We hrst consider the weakly supercritical case (i). Note that for each x, y > 0 

h{x,y) < y^/^). 


Since r_^ and 


are independent, we deduce 


E 


h 



21 


i-v) 


E 


21 


i-v) 






where g{u) ;= E \h «)]. From the inequality of above, we get g{u) < C{u V 

for C > 0 that depends on k, f3 and p. 

Following step by step the proof of Lemma 4.5 in [6], we can deduce that the statement 
also holds for our function g with 5 = 1. Actually in the proof of Lemma 4.5 in [6], the 
authors use the inequality on their statement in order to apply the Dominated Convergence 
Theorem and they split the integral in (4.24) in [6] into two integrals, one over [0,1] and 
another over (1, cx)). In our case, we can take on the integral over [0,1] the function Cu and 
on the integral over (1, cxo) the function and the result will not change. Therefore 


j-3/2 

limf^A > 0) = lim-—-E 

t^QO t^OO LL [Z) 

o roo POO 

= gVr(|„|)I/.(.) I I 

where (j)\ri\ is dehned as in the statement of the Theorem. 

In the remaining two cases we use the following inequalities, which hold by the Mean 
Value Theorem. Let e > 0 then, for each x,y > 0 



^ ^y/y-iy + (^ + l) . (37) 

For the intermediately supercritical case (ii), we note that —rj = 2. From Lemma [T] (with 
p = 1) and Lemma 4.4 in [6] we deduce 


lim \/te^*E 

t —>-00 




On the other hand, from Lemma 4.5 in [6] with g{u) = , we have 


hmf3/2e2iE 

t—)-oo 


1 

(2/®)V/3 


g{u)(p 2 {u)du, 
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where 02 is defined as in the statement of the Theorem. Therefore by the previous limits, 
the independence between r 2 and identity (136|) and inequalities fl37|) . we have that 
for e > 0 the following inequalities hold 


rk\/2 


P‘^a^/nU^{z) 


E 




l/P 1/0-1 
2 


< > 0 ) 

t—)-oo 


< 


z^ky/2 


f3^ay/nU^{z) 


E 




1//3 


(/r 2 + 


Thus our claim holds true by taking limits as e goes to 0. 

Finally, we use similar arguments for the strongly supercritical case (hi). Observe from 
Lemma [T] and the identity in law by Dufresne fl25]) that 


lime-2(^+'^)*E 

t^OO 


1 


[^- 0 + 2 )] ) 


(38) 


where r_(^+ 2 ) is a Gamma r.v. with parameter 
imply 


hmf^/2gr,2t/2^ 
t—>-CX) 


1 

(2/(-^))i//3 


— {t] + 2). If —7] < 2/13, Lemma 4.5 in [ 6 ] 





(39) 


where 4>\r^\ is defined as in the statement of the Theorem. If —rj = 2//3, from Lemma [T] 
and Lemma 4.4 in [ 6 ], we get 


lim 

t^OO 


1 

(24-'?))1//3 


r(i//3) 


(40) 


Next, if —7] > 2/0, from Lemma [Hand the identity in law by Dufresne fl2^ we get 


lim g-2hi//3+D//3E 

t^OO 


1 

(2/GD)i//3 


= E 


pi//3 

-irj+2/0)\ ’ 


(41) 


where r_(^+ 2 // 3 ) is a Gamma r.v. with parameter —(77 + 2/0). Therefore, from the inde¬ 
pendence between r_^ and inequalities fl37)) and the limits in fl38|) . fl39|) . fi40|) and 

dUD, we deduce that for e > 0 , we have 


zk 


l3U,{z 


-E 




E [r_(^+ 2 )] < hme-2(2—> 0) 


< 


t—>-oo 

z^k 

(3UJ^\ 


-E 


,-kzr 


-(^^r_„ + 6 ' 


1 I// 3-1 


E [r_(^+2)]. 


The proof is completed once we take limits as e goes to 0. 


□ 
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5 The immigration case. 

In this section, we introduce continuous state branching processes with immigration in a 
Brownian random environment. In particular, this class of processes is an extension of 
the Cox-Ingersoll-Ross model in a random environment. For simplicity, we introduce such 
class of processes under the assumption that the branching mechanism posses hnite mean. 

Recall that a CB-process with immigration (or CBI-process) is a strong Markov process 
taking values in [0, cxo], where 0 is no longer an absorbing state. It is characterized by its 
branching mechanism. 


= —an + + / (e — 1 + ux)^{dx), 


(0,oo) 


and its immigration mechanism. 


(j){u) = du + / (1 — e “jz/(df). 


where a G M, 7 , d > 0 and 


(x A x^)yu(dx) + / (1 A x)z/(da;) < 00 . 

Jo 


(42) 


It is well-known that if {Yt,t > 0) is a process in this class, then its semi-group is charac¬ 
terized by 

Ex = exp xMi(A) — J 0(Ms)ds|- , for A > 0, 

where Ut solves 

uo{\) = A. 

According to Fu and Li [TU], a CBI-process can be dehned as the unique non-negative 
strong solution of the stochastic differential equation 


Y^ = Yo+ / {d + aY,)ds+ / V^WsdB, 



f-y,- 


zN{ds, dz, dn) -|- 


0 J{0,oo) Jo 



zM{ds, dz), 


0 J{0,oo) 


where B = {Bt,t > 0) is a standard Brownian motion, iV(ds, dz, dn) and M(ds, dz) 
are two independent Poisson random measures with intensities dsyu(dz)dM and dsi/(dz), 
respectively, satisfying the integral condition fH2ll and N is the compensated measure of 
N. The processes B, N and M are mutually independent. 

Motivated from the dehnition of CBBRE, we introduce a continuous state branching 
process with immigration in a Brownian random environment (in short a CBIBRE) as the 
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unique non-negative strong solution of the stochastic differential equation 


'Zit —Zq + I (d -|- aZg) ds I ■\/2'y'^YsdBs -\- 


a 


ZsdB^;'^ 



zN{ds, dz, du) -|- 


0 J (0,oo) J 0 



(43) 


zM(ds, dz), 


0 J{0,oo) 


where a > 0 and {B[^\t > 0) is a standard Brownian motion independent of B and the 
Poisson random measures N and M. Similarly as in the case with no immigration, we 
dehne the auxiliary process 


+ mt, for f > 0, 


where 


a 


m = a 


The following Theorem provides the existence of the CBIBRE as the unique strong solution 
of 


Theorem 5. The stochastic differential equation has a unique non-negative strong 
solution. The process Z = {Zt,t > 0) is a Markov process and its infinitesimal generator 
A satisfies, for every f G 


= ^a^x^f{x) + (xa + d^f'{x) + xj^fix) 


+ / {f{x + z)- f{x))u{dz) +x (/(x + z)- f{x) - zf{x)) fv{dz). 


(44) 


'(0,oo) 


' (0,cxd) 


Furthermore, the process Z, conditioned on iP®, satisfies the branching property and for 
every t > 0 


E, 


exp < —XZte 


-K, 


( 0 ) 


IJ^(O) 


= exp 


-zvt{f),\,K^^'‘) — f (j)(^vt{r, X, K^^'^)e a.s., 


(45) 


where for every (A, 5) G (M+,C(M+)), Vt : s E [0,t] i-G- Vt{s,X,6) is the unique solution of 
the backward differential equation 


A 

ds 


vfis, A, 5) = e^^fioivtis, A, 6)e '^“), vfit, A, 6) = A, 


(46) 


and '0o(A) = '0(A) + «A, for A > 0. 


Proof. We just explain the main steps of the proof, since it follows from similar arguments 
as those used in Theorem 1. The existence of a strong non-negative solution of (145]) 
is analogous (Theorem 2.5 in [8]) and we omitted. The uniqueness implies the strong 
Markov property and by Ito’s formula it is easy to show that the inhnitesimal generator 
of Z satishes fl44|) . 
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The branching property of Zi conditioned on is due to the CBI-process. Let Zt = 


Zte ^ and F G Ito’s formula again guarantees that F{t, Zt) conditioned 


on i^(o) 

is a local martingale if and only if for every f > 0 , 


'0 


t poo 


0 = I { Z.) + dTf (s, Z,) + fe-''--Z.ij,F(s, Z.) | ds 


’ dx'^' 



'0 ^0 
f*t poo 




Zs ( F{s, Zs + ze ) - F{s, Z^) - Zs)ze 


p(d2:)ds 


+ 



0 Jo 


F{s, Zs + ze — F{s, Zs) ) z/(dz)ds. 


( 0 ) 


We choose F(s, x) = exp |—a;nt(s, A, (j){vt{r, X, K^^^)e ■^™)dr| where ^ 4 ( 5 , A, iL®) 

is differentiable with respect to the variable s, non-negative and such that Vt{t, A, = A 

for all A > 0. We observe that F is bounded, therefore, {F{s, Zg), s < t) conditioned on 
RiO) 

is a martingale if and only if 

^vt{s, A, = 7 ^(wt(s, A, ’ 

os 


+ e 


K 


( 0 ) 


^ ^ A, JL® )z ) /i(dz). 


( 0 ) 


which is equivalent to the backward differential equation given in (j46]). The existence and 
unicity of Vt follows from the Picard-Lindelf theorem. 

—ZsVt{s, A, iL®)|, 0 < s < conditioned on K^^'> is a mar- 


Therefore, the process (exp 
tingale, and hence 


E, 


exp 


XZA IJL® 


= exp 


-zvt{0, A, - I 4>{vt{r, A, lA^'>)e \ . 


□ 


We hnish this section with some examples of CBIBRE, all of them related to the stable 
CBBRE. 


Example 4. (Stable CBIBRE) Here we assume that the branching and immigration 
mechanisms are of the form '0(A) = —aA-l-cA^'''^ and 0(A) = kX^, where (3 e (0,1) , c, K > 0 
and a G M. Hence, the stable CBIBRE-process is given as the unique non-negative strong 
solution of the stochastic differential equation 


pt pt pt poo pZs— ^ pt poo 

Zt=ZQ+ / aZgds + a / ZgdB^^^ + / / / ^iV(ds, dz, dw)-|- / / zM{ds,dz), 

Jo Jo Jo Jo Jo Jo Jo 

where is a Brownian motion, and N and M are two independent Poisson random 
measures with intensities 


c/3(/3 + 1 ) 1 

P(1 - (3) z2+/3 


dsdzdw 


and 


k/3 1 

r{l-(3)z^+^ 


dsdz;. 
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Its associated infinitesimal generator A satisfies, for every / G C'^(M+), 


^f(x) = + xaf{x) + f 

Z i (^i - J(o,oo) 


if{x + z) - fix)) 


dz 

^l+/3 


X 


From (H5|) we get the following a.s. identity 

-1//3- 


dz 


(47) 




exp < —XZte 




li^(O) 


= exp l-zi^c I + A"^ 


X 


exp In ^/ScA* e-/”‘-°’ds + 1^ 


(48) 


If we take limits as 2 ; goes to 0, we deduce that the entrance law at 0 of the process 

_44'(0) . . 


(Zje t > 0) satisfies 


Er 


exp < —XZte 




|iF® =exp|-^ln J e'^^^ds + 1 


If we take limits as A goes to 00 in (H8l) . we obtain 

> 0 = 1, for z>0. 


Similarly if we take limits as A goes to 0 in fldHj) . we deduce 

P, (^Zt < 00 = 1, for ^ > 0. 


In other words, the stable CBIBRE is conservative and positive at finite time a.s. 

An interesting question is to study the long-term behaviour of the stable CBIBRE. 
Now, if we take limits as t goes to 00 in (HHj) . we deduce that when m > 0, 


E, 


exp A lim Z^e ' 

L t—^OO 


=E 


exp < 


(3c I + A-^ 


-1//3' 


X exp 


-^In 


(3c 


(3cX^ 




+ 1 


where we recall that /q°° e ^ds has the same law as , with F^ is a Gamma 

r.v. with parameter v. In other words, Zte~^^ ' converges in distribution to a r.v. whose 
Laplace transform is given by the previous identity. 

If m < 0, we deduce 

lim ZfC ‘ = cx), P^ — a.s. 

t^OO 
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We observe that when m = 0, the process oscillates implying that 

lim Zf = oo, P;j — a.s. 

t—^OO 


Example 5. (Stable CBBRE Q-process). Here, we assnme m < —a^. From 
Proposition [3l we dednce that the form of the inhnitesimal generator of the stable CBBRE 
Q-process, here denoted by satishes for / G Dom(£) 


C^f{z) = Cf{z) + z'^a^fiz) 
cf3{(3 -f-1) z 


U'{z) 

U{z) 


r(i-/3) U{z) j. 


f{y + z)- f{z)^ (u{z + y)-U (z)) 


dy 

(,2+/3 ' 


where U was defined as 


U{z) 





kP 



nv - V/g) 

T{y - 2/(3) 


if m = —cr^, 
if m < —cr^. 


Replacing the form of U in the infinitesimal generator in both cases, we get 


+ (a + (^^)zf{z) + 


c(3{P + l) 


+ z 


c(3{(3 -I-1) f 
r(l ~ P) 7(0,00) 


(fix 


r(i-/3) Jo 

-f{z)-xf{z)) 


f{z + 

dx 


X 


fix] 


d^ 


X 


X 


2+P' 


From the form of the inhnitesimal generator of the stable CBIBRE given in flTTD . we deduce 
that the stable CBBRE Q-process is a stable CBIBRE with branching and immigration 
mechanisms given by 


i/j{X) =—{a"^ + a)X + and 0(A) = c(/3-f 1)A^, 


and random environment > 0). 
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